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We  report  one-  and  two-photon  absorption  excitation  energies  and  cross  sections  for  a  series  of 
7-aminocoumarins  using  time-dependent  density  functional  theory  with  various  basis  sets  and 
functionals,  including  exchange-correlation  functionals  using  the  Coulomb-attenuating  method,  to 
evaluate  their  performance  in  the  gas  phase  and  in  solvents.  Except  for  the  results  of  one  functional, 
the  computed  one-photon  excitation  energies  and  transition  dipole  moments  are  in  good  agreement 
with  experiment.  The  range  of  errors  obtained  from  various  functionals  is  discussed  in  detail.  The 
relationship  of  donor  and  acceptor  groups  with  the  one-  and  two-photon  resonances  and  intensities 
is  also  discussed.  ©  2007  American  Institute  of  Physics.  [DOI:  10.1063/1.2464110] 


I.  INTRODUCTION 

Coumarin  derivatives  that  are  formed  by  various  substi¬ 
tutions  of  amino  groups  at  the  7-position  are  widely  used  as 
probes  in  ultrafast  solvation  studies,  as  laser  dyes,1’-  and  as 
nonlinear  optical  chromophores3  in  optoelectronic  devices. 
Thus,  pertinent  optical  properties  for  these  applications  have 
been  extensively  reported  in  the  literature.  Electronic  one- 
photon  absorption  (OPA)  spectra  of  coumarins,  in  particular, 
have  been  characterized  in  high  temperature  vapors,6'7  in  a 
supersonic  jet  expansion,6-8  and  in  a  variety  of  polar  and 
nonpolar  solvents.9-16  Experimental  studies  have  revealed 
the  polar  nature  of  the  first  excited  state  that  has  been  char¬ 
acterized  as  an  intramolecular  charge-transfer  (ICT)  state, 
implying  a  larger  dipole  moment  than  in  the  ground  state. 
The  large  difference  in  dipole  moment  between  the  ground 
and  excited  states  (A/uge)  provides  a  viable  mechanism  for 
the  intense  two-photon  absorption16-19  (TPA)  and  other  non¬ 
linear  processes3'4  that  have  been  observed  and  quantified  for 
a  number  of  coumarin  derivatives. 

Theoretically,  OPA  properties  of  coumarins  in  the  gas 
phase  have  been  reported  in  semiemperical  studies  using 
AMI,20  PM3, 21,22  and  MNDO23  methods.  Recently,  Cave 
et  a/.-4  reported  a  detailed  theoretical  study  of  coumarin  120 
(Cl 20)  and  C151  using  a  wide  range  of  electronic  structure 
methods  to  evaluate  various  ground  and  excited  state  prop¬ 
erties.  They  found  that  using  time-dependent  density  func¬ 
tional  theory-6--8  (TDDFT)  with  hybrid  functionals  gives  ex¬ 
cellent  agreement  with  .S'(r.Sj  experimental  OPA  excitation 
energies.  In  another  survey-  involving  C102,  Cl 52,  Cl 53, 
and  C343,  excellent  agreement  with  experiment  was  also  re¬ 
ported  for  TDDFT  excitation  energies.  Optical  properties  of 
Cl 20, 24  051, 24,30-32  152 A  (C35,  C481), 30,31  and  053 
(Refs.  30  and  31)  have  been  studied  using  TDDFT  with 
continuum,-4  and  explicit30-32  solvent  models.  Recently,  In- 
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grosso  et  al.  reported"  a  solvation  study  of  053  in  super¬ 
critical  fluoroform  using  TDDFT  with  Becke’s 
three-parameter34-36  (B3LYP)  hybrid  functional.  More  re¬ 
cently,  TDDFT  calculations  using  the  conventional 
model-  ’  and  the  state  specific  polarizable  continuum 
model  (PCM)  have  been  carried  out  for  053.  The  successful 
applications  of  TDDFT  with  hybrid  functionals  for  cou¬ 
marins  are  consistent  with  our  previous  studies40  on  the  pre¬ 
diction  of  excitation  energies  involving  charge-transfer  ex¬ 
cited  states.  However,  we  found  that  the  B3LYP 
functional  '4-  6  and  the  hybrid  functional  of  Perdew,  Burke, 
and  Erzerhof4 1 ,4-  (PBEO)  underestimate  the  experimental  ex¬ 
citation  energies  of  a  number  of  ICT  chromophores.  Both 
functionals  were  found  to  overestimate  the  oscillator 
strengths  for  the  first  ICT  excited  states.  The  failures  to  pro¬ 
vide  a  good  description  for  ICT  excited  states  have  been 
attributed  to  the  incorrect  long-range  behavior  of  the  ex¬ 
change  potential.43,44  Thus,  a  detailed  comparison  between 
computed  and  experimental  OPA  and  TPA  excitation  ener¬ 
gies  and  intensities  for  a  series  of  coumarins  was  carried  out 
to  assess  the  accuracy  that  can  be  achieved  with  the  B3LYP, 
PBEO,  and  the  new  long-range-corrected  functional  of 
Handy  and  co-workers.43'44  Previous  TDDFT  studies  on  cou¬ 
marins  focused  on  the  OPA  excitation  energy  but  not  on 
intensity,  and  TPA  spectra  were  not  considered. 

The  aim  of  the  present  work  is  to  investigate  both  the 
OPA  and  TPA  spectra  for  a  series  of  coumarin  derivatives 
with  different  degrees  of  ICT  character  (Figs.  1  and  2)  in  the 
gas  phase  and  in  solvents.  In  7-aminocoumarins,  the  amine 
and  carbonyl  moieties  are  expected  to  act  as  electron  donor 
and  acceptor,  respectively.  Thus,  the  4-methyl  coumarin  de¬ 
rivatives  (Fig.  1)  with  an  increasing  degree  of  alkylation 
measure  the  effects  of  different  donor  groups  on  the  OPA  and 
TPA  spectra.  The  effects  of  the  acceptor  are  also  gauged  by 
4-trifluoromethyl  coumarin  derivatives  (Fig.  2). 

II.  COMPUTATIONAL  METHODS 

The  computation  of  excitation  energies  and  cross  sec¬ 
tions  was  preceded  by  the  determination  of  equilibrium  ge- 
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FIG.  1.  Structure  of  coumarin  120,  1,  480  (102). 

ometries,  and  then  single  point  excited  state  calculations 
were  carried  out  to  obtain  the  vertical  excitation  energies, 
transition  dipoles,  and  oscillator  strengths.  All  structures 
have  been  predicted  using  the  Kohn-Sham43  density  func¬ 
tional  theory  with  the  6-31  IG(d.p)  (Ref.  46)  basis  set  and 
the  B3LYP  hybrid  functional.'4  36  These  structures  were 
verified  to  be  minima  on  the  potential  energy  surface  with 


C152 


CI51 


FIG.  2.  Structure  of  coumarin  152,  152A,  522,  153,  307,  and  151. 


harmonic  frequency  calculations  using  the  GAUSSIAN  03 
program.47  TDDFT  calculations  were  carried  out  using  the 
6-31  lG(d,p),  6-311  +  +G (d,p),  and  aug-cc-pVTZ  (for  se¬ 
lected  cases)  basis  sets  with  the  B3LYP,  PBE0,  and  the  new 

as  implemented  in  the 
functional,  called 
CAMB3LYP,  combines  properties  of  B3LYP  and  the  long- 
range  correction  for  the  exchange  potential  by  an  Ewald  split 


hybrid  functional  of  Yanai  et  a/.,43 

48 

DALTON  program.  This  hybrid 


of  the  r ,  2  operator  into43 


r-1 . 
'  12  ■ 


1  -  [a  +  /3  erf (/rr12)]  a  +  f3  erf(/rr12) 


r\2 


r  12 


(i) 


where  p,  a,  and  (i  are  empirical  parameters.  Yanai  et  al. 
found  that  the  CAMB3LYP  functional  yields  the  smallest 
mean  absolute  error  (MAE=2.53  kcal/mol)  with  /x= 0.33, 
q,=0.19,  and  /3=0.46  (a+/3=0.65)  in  atomization  energies 
for  the  G2  set.  This  is  slightly  better  than  the  corresponding 
MAE  of  2.68  kcal/mol  for  B3LYP,  which  is  nearly  equiva¬ 
lent  to  setting  a= 0.2  and  /3=  0.0  in  CAMB3LYP.44  However, 
the  ICT  excitation  energies  for  a  model  dipeptide  are  greatly 
improved,  when  compared  to  the  CASPT2  results.  In  addi¬ 
tion,  we  tested  a  modified  CAMB3LYP  (mCAMB3LYP) 
functional  with  a  two-parameter  Ewald  split  of  the  rj2  opera¬ 
tor  as 


-l 

r12  : 


1  -  [a+  a  erf(/rr 12)]  a  +  a  erf(/rrp) 

- + - — . 


G  2 


r  12 


(2) 


where  p.  and  a  are  kept  unchanged. 

To  account  for  the  solvent  effects,  TDDFT  calculations 
were  carried  out  with  the  nonequilibrium  PCM  (Refs.  37  and 
38)  and  the  self-consistent  reaction  field  (SCRF)  model  of 
Mikkelsen  et  al.,  which  uses  a  spherical  solute  cavity 
(SCRF-S).  The  PCM  solute  cavities  were  carried  out  using 
the  united  atom  topological  model  and  empirical  solvent  ra¬ 
dii,  as  implemented  in  the  GAUSSIAN  03  program  47  The  radii 
(a0)  f°r  the  SCRF-S  were  obtained  using  the  largest  inter¬ 
atomic  distances  of  the  solutes  with  added  adjustments  for 
van  der  Waals  interactions.  The  SCRF-S  calculations  were 

48 

carried  out  using  the  DALTON  program. 

The  TPA  cross  sections  can  be  obtained  by  relating  the 
absorption  rate40'50  to  the  TPA  transition  probability,  which 
was  first  derived  by  Goppert-Mayer  (GM)  using  second  or¬ 
der  perturbation  theory.  1  The  two-photon  absorptivity  <5/o 
between  the  ground  state  (0)  and  the  final  state  /  is  written 


as 


40,52,53 


8  77 


<5/o(£1’£2)  =  gjQ(El,E2)\Sf0{El,E2)\2,  (3) 


/ 


where  c  is  the  speed  of  light,  h  is  Planck’s  constant,  Ek  the 
photon  energies,  and  g/0(£j,£2 )  is  the  normalized  line  shape 
function  for  the  ground  to  /  transition,  which  can  be  repre¬ 
sented  by  a  Gaussian  or  a  Lorentzian.  The  two-photon  matrix 
elements  Sj-0(E1,E2)  were  obtained  from  single  (SRQR)  or 
double  (DRQR)  residue  of  the  quadratic  response  function  of 
the  electric  dipole  operator34  within  the  framework  of 
TDDFT.53  TPA  cross  section  calculations  with  DRQR  are 

carried  out  using  the  sum-over-state  (SOS)  method.  These 

48 

calculations  were  carried  out  with  the  DALTON  program. 
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The  excited  state  dipoles  are  computed  using  the  DRQR 
method  or  the  finite  field  method  with  an  applied  field  of 
±0.0001  a.u. 

III.  RESULTS  AND  DISCUSSION 

The  OPA  results  are  briefly  summarized  in  Table  1  while 
all  the  results  are  organized  in  an  expanded  version  of  the 
table  in  the  EPAPS  (Ref.  68)  Supporting  Information  Sec¬ 
tion.  The  TPA  results  are  tabulated  in  Tables  II-IV.  Table  I 
lists  the  computed  excitation  energies,  transition  dipole  mo¬ 
ments  (/u.ge),  and  oscillator  strengths  if)  along  with  available 
experimental  values.  The  experimental  energies  are  obtained 
from  the  maxima  of  the  absorption  band,  which  were  inte¬ 
grated  (using  Gaussian  line  shapes)  to  obtain  the  transition 
moments  and  oscillator  strengths  in  cases  where  the  absolute 
intensities  were  available.  We  note  that  the  gas-phase  B3LYP 
excitation  energies  for  Cl 20  and  C151  and  the  corresponding 
PBE0  values  for  C152  and  C153  are  identical  to  those  re¬ 
ported  by  Castner  and  co-workers.24"9  Table  I  also  lists  the 
OPA  spectral  properties  in  ethanol/methanol  and  cyclohex¬ 
ane.  These  solvents  are  chosen  due  to  their  widespread  use  as 
representatives  for  polar  and  nonpolar  solvents  for  cou- 
marins.  We  begin  with  the  OPA  spectra,  and  then  discuss  the 
TPA  results,  starting  with  the  gas-phase  spectra,  the  influence 
of  basis  sets,  theoretical  methods,  and  the  effects  of  solvents 
and  alkylation.  Analysis  of  the  basis  set  effect  is  limited  to  a 
few  systems,  which  are  found  to  be  small.  Thus,  the  results 
of  the  6-31 1G (d,p)  basis  set  are  used  for  discussion  through¬ 
out,  unless  noted  otherwise. 

A.  One-photon  absorption 

We  first  compare  the  gas-phase  results  that  are  available 
for  C152,  C152A,  C522,  and  C153  (Table  I).  For  these  cou¬ 
marin  derivatives,  TDDFT  excited  state  calculations  are  car¬ 
ried  out  at  the  B3LYP/6-3I  lG(d,p)  geometry  using  both 
the  6-311G(<7,/?)  and  6-311  +  +G(c/,/7)  basis  sets  to  gauge 
the  effects  of  diffuse  functions.  For  CAMB3LYP  (using  the 
recommended  a=0.19  and  /j=0.46  parameters)  and  the  two 
widely  used  B3LYP  and  PBE0  functionals,  the  TDDFT  ex¬ 
citation  energies  are  slightly  affected  by  the  basis  set.  The 
augmentation  of  diffuse  functions  on  hydrogen  and  heavy 
atoms  lowers  the  excitation  energies  by  about  0. 1  eV  for  the 
,8'0-.S’ ]  ICT  energies,  as  predicted  by  B3LYP  and  PBE0.  How¬ 
ever,  the  corresponding  changes  are  found  to  be  quite  small 
(0.02  eV)  for  the  CAMB3LYP  functional.  Overall,  the  ICT 
experimental  excitation  energies  are  well  reproduced  with 
the  6-31 1G {d,p)  basis  set,  which  yields  the  MAEs  of  0.03 
and  0.08  eV  for  PBE0  and  B3LYP,  respectively.  The  corre¬ 
sponding  MAEs  of  0.16  and  0.05  eV  are  found  for  the 
6-31 1  +  +G(d,p)  basis  set,  which  has  negligible  differences 
in  Sq-Si  excitation  energies  from  the  larger  aug-cc-pVTZ  ba¬ 
sis  set  (see  Table  I).  Note  that  B3LYP  consistently  underes¬ 
timates  the  experimental  excitation  energies.  The  predicted 
CAMB3LYP/6-31  \G{d,p)  ICT  excitation  energies  yield  a 
MAE  of  0.32  eV.  This  is  much  larger  than  the  corresponding 
MAE  obtained  with  the  PBE0  and  B3LYP  functionals.  As 
illustrated  for  C152,  the  large  deviations  between  B3LYP  and 
CAMB3LYP  can  be  attributed  to  the  value  of  f3  [see  Table  IS 


in  the  EPAPS  (Ref.  68)  Supporting  Information  Section], 
since  /3=  0  for  B3LYP  (a=0.19  nearly  the  same  as  in 
B3LYP).  Thus,  the  a  and  (i  parameters  must  be  reexamined 
for  both  the  ground  (e.g.,  the  G2  set)  and  excited  states. 
However,  the  effects  of  f3  on  the  ground  state  energetics  are 
found  to  be  small.  The  MAE  of  atomization  in  the  G2  set  for 
f3=  0  and  /3=0.46  changes  from  2.68  to  2.53  kcal/mol, 
respectively.43  Thus,  the  corresponding  error  for  /3=a= 0.19 
[Eq.  (2)]  is  likely  within  this  range.  The  small  compromise  in 
the  ground  state  energetics  results  in  a  significant  improve¬ 
ment  in  the  predicted  S0-S ,  ICT  energies,  which  have  MAEs 
of  0.10  and  0.06  eV  (without  C153  due  to  convergence  prob¬ 
lem)  for  the  6-311G (d,p)  and  6-31 1  +  +G(d,p)  basis  sets, 
respectively.  It  is  worth  noting  that  excitation  energies  com¬ 
puted  with  the  nonhybrid  BLYP  functional  and  Tamm- 
Dancoff  approximation  are  computationally  more  economi¬ 
cal  but  do  not  offer  comparable  accuracy  delivered  by 
TDDFT  with  the  hybrid  B3LYP,  PBE0,  or  mCAMB3LYP 
functional.  The  predicted’0  '1  gas-phase  S0-Sl  energies  of 
2.99  and  2.90  eV  for  C152A  (C35)  and  C153,  respectively, 
are  significantly  lower  than  the  corresponding  experimental 
values6  of  3.56  and  3.37  eV.  For  051,  the  BLYP  and  SAOP 
gas-phase  excitation  energies  of  3.32  (Refs.  30  and  31)  and 
3.24  eV,  respectively,  were  reported.  These  predicted  val¬ 
ues  are  much  lower  than  the  experimental  value  of  3.70  eV 
in  n-hexane. 1  2 

We  now  discuss  the  trends  observed  with  the  predicted 
results  in  the  gas  phase,  which  do  not  appear  to  be  dependent 
upon  the  choice  of  basis  and  functionals.  For  example,  the 
B3LYP  and  PBE0  (0.17  eV)  Sq-Sj  shifts  for  C152-C153  are 
identical  for  the  6-31 1G**  and  6-31 1  +  +G**  basis  sets.  The 
predicted  value  of  0.17  eV  is  in  good  agreement  with  corre¬ 
sponding  experimental  value  of  0.22  eV.  All  four  functionals 
predict  the  Sq-S\  energy  to  shift  to  the  red  with  an  increasing 
degree  of  alkylation,  C120<C1  <C480.  A  trifluorine  substi¬ 
tution  at  the  4-methyl  group  gives  pronounced  redshifts  that 
increase  in  magnitude  with  a  higher  degree  of  alkylation 
(C151  <C307<C152<  152A<C522<C153),  in  agree¬ 
ment  with  the  experimental  data6  for  the  latter  four  cou- 
marins.  We  should  note  that  the  absorption  maxima  of  Cl 52 
and  C152A  listed  in  Table  I  of  Ref.  6  must  be  interchanged 
for  them  to  be  consistent  with  their  corresponding  absorption 
spectra  given  in  Fig.  3  of  their  paper. 

The  continuum  solvation  models  do  not  treat  specific 
solvent-solute  interactions  such  H  bonding,  which  is  impor¬ 
tant  in  water8"4  and  methanol50  but  not  in  a  cyclohexane 
solution.  For  the  6-311G (d,p)  basis  set,  the  SCRF-S  model 
of  cyclohexane  overestimates  the  experimental  S0-Sl  ener¬ 
gies  by  0.05,  0.22,  and  0.42  eV  for  the  B3LYP, 
mCAMB3LYP,  and  CAMBLYP,  respectively.  Using  the 
same  basis  set,  the  MAEs  for  PCM  with  PBE0  and  B3LYP 
are  0.03  and  0.11  eV,  respectively.  These  excitation  energies 
do  not  appear  to  be  strongly  affected  by  the  addition  of  dif¬ 
fuse  functions,  which  are  found  to  lower  the  energies  by 
0.05-0.07  eV  for  C120  (see  Table  I).  Thus,  the  6-311G (d,p) 
basis  set  seems  adequate  for  predicting  S^-Sy  energies.  Note 
that  the  errors  for  B3LYP  are  slightly  improved  compared  to 
those  obtained  in  the  gas  phase  due  to  the  effects  of  compen¬ 
sation  in  excitation  energy  by  the  SCRF-S  method  and  the 
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TABLE  I.  6-31  \G(d,p)  excitation  energy  ( E ,  eV),  oscillator  strength  (/),  and  transition  dipole  moment  (/x01,D)  for  coumarin  derivatives  compared  with 
experiment.  The  values  in  parentheses  and  brackets  are  computed  with  the  6-31 1  +  +G(d,p)  and  the  aug-cc-pVTZ  basis  sets,  respectively. 


Theor. 

Expt. 

System/method 

E(SJ 

/ 

Mm 

Energy  (/t0i)  [/] 

Solvent 

Cl  20 

CAMB3LYP 

4.29  (4.24) 

0.373  (0.368) 

4.79  (4.79) 

PBE0 

4.09  (4.04)  [4.03] 

0.348  (0.363)  [0.355] 

4.74  (4.87)  [4.82] 

B3LYP 

4.00  (3.94)  [3.93] 

0.329  (0.344)  [0.338] 

4.66  (4.80)  [4.76] 

EtOH 

SCRF-S-CAMB3LYPa 

4.07  (4.00) 

0.479  (0.497) 

5.57  (5.72) 

SCRF-S-mCAMB3LYPa 

3.90  (3.83) 

0.433  (0.451) 

5.41  (5.57) 

PCM-PBE0 

3.81  (3.74) 

0.456  (0.477) 

5.61  (5.80) 

SCRF-S-B3LYPa 

3.76  (3.69) 

0.390  (0.406) 

5.23  (5.39) 

PCM-B3LYP 

3.71  (3.64) 

0.424  (0.445) 

5.49  (5.68) 

3.5 lb  (5.54)  [0.3 1 8]cd 

Water 

SCRF-S-CAMB3LYPe 

(3.99) 

(0.495) 

(5.72) 

SCRF-S-mCAMB3LYPe 

(3.81) 

(0.447) 

(5.56) 

SCRF-S-B3LYPe 

(3.68) 

(0.401) 

(5.36) 

PCM-B3LYP 

3.70  (3.63) 

0.418  (0.439) 

5.46  (5.65) 

3.65  (4.85)  [0.325]f 

Cyclohexane 

SCRF-S-CAMB3LYPe 

4.19  (4.13) 

0.467  (0.477) 

5.42  (5.51) 

SCRF-S-mCAMB3LYPe 

4.02  (3.97) 

0.434  (0.446) 

5.33  (5.44) 

SCRF-S-B3LYPg 

3.89  (3.84) 

0.401  (0.415) 

5.22  (5.33) 

PCM-PBE0 

3.93  (3.87) 

0.472  (0.495) 

5.63  (5.81) 

PCM-B3LYP 

3.84  (3.77) 

0.443  (0.468) 

Cl 

5.52  (5.72) 

3.76b 

B3LYP 

3.72  (3.66)  [3.65] 

0.401  (0.409)  [0.405] 

5.33  (5.43)  [5.41] 

PCM-B3LYP 

3.47 

0.490 

6.10 

3.32b  (5.63)  [0.399]dh 

EtOH 

PCM-B3LYP 

3.57 

0.510 

C480 

6.13 

3.54f 

Cyclohexane 

B3LYP  ( syn ) 

3.60 

0.347 

5.04 

PCM-B3LYP  (anti) 

3.33 

0.447 

5.95 

MeOH 

PCM-B3LYP  (syn) 

3.34 

0.433 

5.85 

3.19  (5.84)  [0.412]d>1 

MeOH 

PCM-B3LYP 

3.45 

0.455 

C151 

5.89 

3.43J 

Cyclohexane 

B3LYP 

3.73  (3.66)  [3.67] 

0.312  (0.320)  [0.317] 

4.70  (4.80)  [4.77] 

TDA-BLYP*1 

3.32 

SAOP/TZP1 

3.24 

PCM-B3LYP 

3.38 

0.393 

5.54 

3.24  (4.97)  [0.304]c 

EtOH 

MD-QMMM-SAP/TZP1 

3.07 

3.70m 

n-hexane 

PCM-B3LYP 

3.53 

0.411 

C307 

5.54 

3.56° 

Cyclohexane 

B3LYP 

3.55 

0.352 

5.12 

PCM-B3LYP 

3.27 

0.431 

5.89 

3.15  (5.52)  [0.364]° 

MeOH 

PCM-B3LYP 

3.38 

0.451 

C152 

5.93 

3.37J 

Cyclohexane 

B3LYP 

3.48  (3.40) 

0.352  (0.348) 

5.17  (5.20) 

3.59p 

PCM-B3LYP 

3.20 

0.427 

5.93 

3.13f  (5.31)  [0.334]q 

EtOH 

C152A  (C35) 

B3LYP 

3.44  (3.35) 

0.375  (0.372) 

5.37  (5.41) 

3.56p 

TDA-BLYPk 

2.99 

PCM-B3LYP 

3.16 

0.457 

6.18 

3.10)  (5.88)  [0.408]d'm 

EtOH 

PCM-B3LYP 

3.27 

0.473 

C522 

6.17 

3.30> 

Cyclohexane 

B3LYP 

3.42  (3.34) 

0.370  (0.368) 

5.34  (5.39) 

3.45p 

PCM-B3LYP 

3.14 

0.463 

6.24 

3.05J 

EtOH 

PCM-B3LYP 

3.24 

0.475 

6.21 

3.25J 

Cyclohexane 

4 
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TABLE  I.  ( Continued .) 


System/method 

£(St) 

Theor. 

/ 

Aha 

Expt. 

Energy  (yu.0t)  [f] 

Solvent 

TDA-BLYPk 

B3LYP 

PCM-B3LYP 

PCM-B3LYP 

2.90 

3.31  (3.23) 

3.04 

3.15 

0.331  (0.333) 
0.422 

0.434 

C153 

5.14  (5.21) 

6.05 

6.03 

3.37p  (5.51)  [0.37]r 

2.95'  (5.71)  [0.373]” 
3.15'  (5.84)  [0.408]5 

EtOH 

Cyclohexane 

aa0=  10.00  a.u.,  £  =  24.55,  n2=  1.847. 
bReference  14. 

“Reference  16. 

^Reference  61. 

ea0=10.0  a.u.,  £=78.39,  and  n2=  1.776. 
'Reference  18. 

8a0=11.00  a.u..  £=24.55,  andn2=1.847. 
hReference  65. 

‘Reference  19. 

'Reference  9. 


kReference  31. 

'Reference  32,  excitation  energy  in  n-hexane  was  obtained  from  an 
average  of  400  snapshots  of  the  molecular  dynamics  (MD)  simulations. 
“Reference  12. 

“Reference  64. 

“Reference  69. 

pExperimental  gas-phase  .Vq-.S).  Reference  6. 
qReference  60. 

“Reference  7. 

“Reference  11. 


B3LYP  functional.  This  might  be  attributed  to  the  tendency 
to  overestimate  and  underestimate  the  excitation  energy  by 
the  SCRF-S  model  and  the  B3LYP  functional,  respectively. 
Consequently,  the  mCAMB3LYP  and  CAMBLYP  energies 
that  are  overestimated  in  the  gas  phase  are  exacerbated  in  the 
SCRF-S  model.  For  B3LYP,  the  SCRF-S  excitation  energies 
are  larger  (by  0.03-0.07  eV)  compared  to  the  corresponding 
PCM  (with  MAE  of  0.03  eV)  values.  This  suggests  that  ex¬ 
citation  energies  of  mCAMB3LYP  and  CAMBLYP  are 
likely  improved  with  the  PCM  and  the  use  of  diffuse  func¬ 
tions,  assuming  the  reduction  in  excitation  energy  observed 
in  C120  holds  for  other  coumarins.  Thus,  using  PCM  with 
more  flexible  basis  sets  might  bring  the  mCAMB3LYP  er¬ 
rors  closer  to  those  obtained  with  PBE0  (MAE  of  0.11  eV). 
However,  it  is  unlikely  that  using  PCM  with  larger  basis  sets 
would  bring  the  CAMB3LYP  errors  to  an  acceptable  range 
for  these  systems.  Overall,  we  found  that  the  predicted 


B3LYP  and  PBE0  S’q-xS' x  energies  in  cyclohexane  are  in  ex¬ 
cellent  agreement  with  experiment.  The  errors  for 
mCAMB3LYP  and  CAMB3LYP  are  slightly  larger  than 
those  obtained  in  the  gas  phase.  The  trends  predicted  for  the 
gas  phase  S0-S1  energy  to  shift  to  the  red  with  an  increasing 
degree  of  alkylation  for  the  nonfluorinated  (C120<C1 
<C480)  and  for  the  trifluorinated  series  (C151<C307 
<C152<  152A<C522<C153)  also  hold  true  in  cyclohex¬ 
ane,  and  are  confirmed  by  a  collection  of  experimental  val¬ 
ues  in  different  studies. 

For  ethanol,  the  SCRF-S  model  and  the  PCM  appear  to 
overestimate  S0-Sl  energies,  suggesting  the  H-bonding  ef¬ 
fects  are  significant,  especially  in  Cl 20  and  C151.  The  over¬ 
all  errors  for  PCM,  however,  remain  low  for  B3LYP 
(0.11  eV)  and  PBE0  (0.21  eV).  The  MAEs  using  the 
SCRF-S  model  are  0.13,  0.27,  and  0.48  eV  for  the  B3LYP, 


TABLE  II.  TPA  cross  section  (in  GM,  using  FWHM=0.2  eV)  and  excited  state  maximum  (in  parentheses) 
for  C120.  Experimental  TPA  cross  section  of  19  GM  (not  a  maximum  value)  measured  at  3.29  eV  in  ethanol, 
(Ref.  16);  experimental  TPA  cross  section  of  14  GM  (not  a  maximum  value)  measured  at  4. 13  eV  in  water  (Ref. 
18). 


Method 

SRQR 

DRQR-SOS  (5) 

A 

So 

A 

CAMB3LYP/6-31G* 

5.1(4.34) 

3.6(4.34) 

6.8(4.35) 

4.7(4.35) 

PBEO/6-31G* 

6.4(4.13) 

4.4(4.13) 

6.9(4.14) 

4.7(4.14) 

B3LYP/6-31G* 

6.6(4.03) 

4.6(4.04) 

6.9(4.05) 

4.8(4.05) 

CAMB3LYP/6-311G** 

5.2(4.30) 

3.7(4.30) 

6.5(4.31) 

4.4(4.31) 

PBEO/6-311G* 

6.5(4.10) 

4.5(4.10) 

6.5(4.12) 

4.5(4.11) 

B3LYP/6-311G** 

6.7(4.01) 

4.7(4.01) 

6.8(4.02) 

4.7(4.01) 

CAMB3LYP/6-3 1 1  +  +  G” 

5.5(4.25) 

3.9(4.25) 

5.9(4.26) 

4.0(4.26) 

PBEO/6-311  +  +G** 

6.9(4.06) 

4.8(4.06) 

6.8(4.07) 

4.7(4.07) 

B3LYP/6-311++G” 

7.1(3.96) 

5.0(3.96) 

7.0(3.97) 

4.9(3.97) 

SCRF-S-B3LYP/6-3 1 1  +  +  G**a 

13.9(3.69) 

9.5(3.69) 

17.4(3.70) 

11.7(3.69) 

SCRF-S-CAMB3LYP/6-3 1 1  +  +  G”a 

13.2(4.00) 

9.0(4.00) 

15.5(4.01) 

10.6(4.01) 

SCRF-S-B3LYP/6-311  +  +  G**b 

14.1(3.70) 

9.6(3.70) 

17.6(3.72) 

11.9(3.71) 

SCRF-S-CAMB3LYP/6-3 1 1  +  +  G"b 

13.8(4.01) 

8.9(4.01) 

15.3(4.03) 

10.5(4.02) 

aa0=  10.00  a.u.,  £=78.39,  and  n2=  1.776  (water). 
b«o=10  a.u.  £  =  24.55,  and  «2=  1.847  (EtOH) 
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TABLE  III.  Ground  and  excited  state  dipole  moments  (in  D)  and  TPA  cross  sections  for  C120  using  the 
two-state  approximation. 


Method 

Mo 

Msi 

Ms2 

Ms3 

Solvent 

CAMB3LYP/6-31G* 

6.3 

8.3 

7.5 

2.6 

4.8 

PBEO/6-31G* 

6.4 

9.0 

8.3 

2.4 

7.5 

B3LYP/6-31G* 

6.3 

9.1 

8.4 

2.2 

8.2 

CAMB3LYP/6-311G* 

6.4 

8.3 

7.4 

2.6 

5.0 

PBEO/6-311G* 

6.4 

8.9 

8.1 

2.3 

7.6 

B3LYP/6-311G** 

6.4 

9.1 

8.2 

2.2 

8.3 

CAMB3LYP/6-3 1 1  +  +  G” 

6.7 

8.5 

7.4 

2.6 

4.7 

PBEO/6-311  +  +G** 

6.6 

9.0 

8.3 

2.2 

7.3 

B3LYP/6-311  +  +  G” 

6.7 

9.2 

8.3 

2.1 

8.0 

Expt.c 

7.9 

EtOH 

PCM-PBE0/6-31G* 

9.4 

14.5 

12.3 

3.7 

33.6 

PCM-B3LYP/6-31G* 

9.4 

14.5 

12.4 

3.7 

33.0 

PCM-PBE0/6-311G** 

9.6 

14.7 

12.0 

11.7 

34.9 

PCM-B3LYP/6-311G" 

9.6 

14.8 

12.2 

11.4 

34.4 

SCRF-S-CAMB3LYP/6-3 1 1  +  +  G**a 

9.5 

12.6 

10.5 

11.4 

14.5 

SCRF-S-B3LYP/6-3 1 1  +  +  G**b 

9.6 

13.4 

11.2 

10.5 

19.1 

Water 

PCM-B3LYP/6-311G" 

9.8 

15.2 

12.4 

11.7 

36.5 

SCRF-S-CAMB3LYP/6-3 1 1  +  +  G**a 

9.7 

12.8 

10.7 

11.7 

14.6 

PCM-B3LYP/6-3 1 1  +  +  G** 

10.4 

15.9 

13.0 

12.3 

42.0 

SCRF-S-B3LYP/6-3 1 1  +  +  G**a 

9.8 

13.6 

11.4 

10.8 

18.9 

aa0=  10.00  a.u.  e=78.39,  and  n2=  1.776  (water). 
ba0=  10.00  a.u.  e=24.55,  and  n2=  1.847  (etOH). 

Estimated  from  the  Lippert-Mataga  relationship  using  fj.g  value  of  6.0  D  and  r  value  of  3.41  A  (Ref.  14). 


mCAMB3LYP,  and  CAMBLYP,  respectively.  The  latter  two 
functionals  overestimate  all  the  .S’q-.S  [  excitation  energies, 
with  Cl 20  having  the  largest  error.  In  fact,  the  continuum 
solvation  models  fail  to  predict  the  observed  blueshift  in  the 
excitation  energy  in  water  with  respect  to  EtOH  (see  Table  I). 
This  might  be  attributed  to  the  H  bonding,  which  is  not  ac¬ 
counted  for  in  the  continuum  models  of  solvation.  In  Cl 20 
(and  051),  water  and  alcoholic  solvents  can  form  hydrogen 
bonds  with  the  nitrogen  (A)  or  oxygen  (B)  at  the  amine  and 
the  carbonyl  groups,  respectively,  and/or  accepting  a  H  bond 
at  the  amine  hydrogens  (C).56  5K  The  B  and  C  interactions 
strongly  stabilize  the  positive  and  negative  charges  on  the 
amino  nitrogen  and  carbonyl  oxygen,  respectively,  in  ICT 
excited  state  of  020,  leading  to  a  redshift  in  the  ICT  exci¬ 
tation  energy  in  the  more  acidic  water  compared  to  alcohols. 
The  A  and  B  interactions  are  present  in  all 
7-aminocoumarins,  while  the  C  interactions  are  unique  to  the 
nonalkylated  and  the  /V-alkylated  (not  /V-dialkylated)  sys¬ 
tems.  However,  the  blueshift  in  the  Sq-S\  excitation  energy  in 
water  with  respect  to  alcohols  has  been  attributed  to  the  A 
interaction  that  better  stabilizes  the  ground  state  than  the  ICT 
excited  state.  ’  Interestingly,  previous  DFT  studies  of 
051  with  one  water  molecule  did  not  locate  a  stable  A-type 
conformer  in  the  ground  state.  Therefore,  calculations  of 
020  in  water  and  methanol  clusters  embedded  in  a  dielec¬ 
tric  continuum  are  important  to  shed  more  light  on  the  effects 
of  H  bonding  on  the  solution  spectra.  Thus,  we  have  carried 
out  PCM-TDDFT  calculations  on  020  with  a  small  number 
of  water  and  methanol  molecules.  Although  some  conform- 
ers  with  the  A  and  C  interactions  possess  blueshifts  in  exci¬ 
tation  energies,  the  computed  Boltzmann  averaged  excitation 


energies  (neglecting  entropy)  are  similar  for  water  and 
methanol  [see  Table  2S  in  the  EPAPS  (Ref.  68)  Supporting 
Information  Section].  The  inclusions  of  entropy,  which  re¬ 
quires  more  sophisticated  treatments  for  low  frequency 
modes  beyond  the  so  called  rigid  rotor  and  harmonic  oscil¬ 
lator  approximation,  are  under  consideration.  Calculations 
for  larger  clusters  using  Monte  Carlo  optimization59  are  also 
being  carried  out  and  will  be  reported  in  the  future. 

We  now  consider  the  OPA  intensity,  which  can  be  mea¬ 
sured  with  the  fjLge  or  /  values  that  are  available  for  Cl 52 
(Ref.  60)  and  C153  (Ref.  11)  in  a  number  of  solvents.  For 
C153,  the  emission  and  absorption  transition  moments  were 
reported  to  be  the  same  to  within  the  experimental 
uncertainty.11  For  other  systems,  the  /xgc  and  /  values  in  so¬ 
lutions  were  obtained  by  using  published  maxima  of  absorp¬ 
tion  coefficients  and  corresponding  full  width  at  half  maxima 
(FWHM),  and  by  integrating  the  appropriate  spectral  bands, 
which  were  approximated  by  Gaussians.  Although  gas-phase 
absorption  spectra  for  a  number  of  coumarins  have  been  re¬ 
ported,  the  corresponding  absolute  intensities  were  not  mea- 
sured,  except  for  the  /  value  of  0.37  for  Cl 53,  which  is 
within  the  range  (0.33-0.44)  of  the  computed  values.  For 
other  systems,  our  comparison  between  theory  and  experi¬ 
ment  is  largely  limited  to  values  obtained  in  alcoholic  solu¬ 
tions.  Experimental  /  values  in  cyclohexane  and  water  are 
also  reported  for  some  systems. 

The  value  of  /  for  C120  in  water  (0.325)  was  obtained 
by  digitization  of  the  OPA  spectrum  reported  by  Overway 
and  Lytle. lis  The  corresponding  value  in  ethanol  (0.318)  was 
obtained  from  the  spectral  data  reported  by  Fisher  et  a/.16 
and  Reynolds  and  Drexhage.61  These  experimental  oscillator 
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TABLE  IV.  Computed  maximum  TPA  cross  section  (in  GM,  using  FWHM=0.2  eV)  and  corresponding  tran¬ 
sition  energy  (in  parentheses)  compared  with  the  experiment.  Single  denotes  “not  a  maximum  value.” 


System/method 

SQR 

Expt. 

A 

S 

SG  (TS)a 

Cl  20 

B3LYP/6-311  +  +  G" 

7.1  (3.96) 

5.0  (3.96) 

SCRF-S-mCAMB3LYP/6-31 1  +  +  G**b 

13.7  (3.82) 

9.3  (3.82) 

SCRF-S-CAMB3LYP/6-31 1  +  +G"b 

13.2  (4.00) 

9.0  (4.00) 

SCRF-S-B3LYP/6-3 1 1  +  +  G**b 

14.1  (3.70) 

9.6  (3.70) 

14  (4.13)°’’ 

4a-d 

SCRF-S-CAMB3LYP/6-3 1 1  +  +  G**e 

13.8  (4.01) 

8.9  (4.01) 

SCRF-S-mCAMB3LYP/6-3 1 1  +  +  G**e 

13.8  (3.84) 

9.4  (3.84) 

SCRF-S-B3LYP/6-3 1 1  +  +  G**e 

13.9  (3.69) 

9.5  (3.69) 

19  (3.29)f’* 

Cl 

B3LYP/6-311G** 

18.6  (3.73) 

12.7  (3.73) 

SCRF-S-CAMB3LYP/6-311G**® 

26.2  (3.91) 

17.8  (3.91) 

SCRF-S-mCAMB3LYP/6-311G**g 

28.0  (3.70) 

19.2  (3.70) 

SCRF-S-B3LYP/6-31  lG**g 

28.6  (3.53) 

19.6  (3.53) 

104  (3.29)f,‘ 

8a,h 

C480  (Cl 02) 

B3LYP/6-311G**  ( syn ) 

14.9  (3.61) 

10.2  (3.61) 

SCRF-S-CAMB3LYP/6-311G”  (syn)1 

21.7  (3.79) 

14.9  (3.79) 

SCRF-S-mCAMB3LYP/6-31 1G”  (syiif 

22.1  (3.58) 

15.2  (3.58) 

SCRF-S-B3LYP/6-311G**  (syn)1 

21.7  (3.42) 

14.8  (3.42) 

160  (3.17y 

9ak 

C151 

B3LYP/6-311G” 

10.3  (3.74) 

7.1  (3.74) 

SCRF-S-B3LYP/6-31  lG**g 

14.0  (3.45) 

9.5  (3.45) 

47  (3.28)f’* 

14a-‘ 

SCRF-S-mCAMB3LYP/6-311G”g 

15.4  (3.60) 

10.5  (3.60) 

SCRF-S-CAMB3LYP/6-311G**® 

16.8  (3.79) 

11.5  (3.79) 

C307 

B3LYP/6-311G** 

15.4  (3.56) 

10.5  (3.56) 

SCRF-S-CAMB3LYP/6-311G**1 

29.0  (3.55) 

19.8  (3.55) 

SCRF-S-mCAMBSLYP/e-SllG**1 

27.1  (3.34) 

18.5  (3.34) 

SCRF-S-B3LYP/6-31 1G*“ 

25.1  (3.17) 

17.1  (3.17) 

19 77  (3.20)” 

aTPA  cross  sections  are  obtained  with  experimental  /a,t|.  in  Table  I.  A/zge  from  time-resolved  microwave  dielec¬ 
tric  absorption  measurements  in  benzene  (Ref.  62),  computed  angles  (C120=160°,  0  =  165°,  C480=164°,  and 
051  =  164°)  between  yu.ge  and  A and  Gaussian  line  shapes  (FWHM=0.2  eV)  using  the  two-state  approxi¬ 
mation. 

b«0=10  a.u.,  £=78.39,  and  «2=  1.776  (water). 

“Experimental  TPA  cross  section  measured  in  water  (Ref.  18). 

°A/xg,.=2.4  I).  using  A/zg(,=4. 1  -6.2  I)  (electro-optic)  (Ref.  62)  gives  <5c(TS)  =  ll-26  GM;  and  using  A yu,ge 
=  1.90  D  (Ref.  14)  (solvatochromism)  gives  5C(TS)  =  2  GM;  and  using  A/ige=3.97  D  (Ref.  66)  (solvato- 
chromism)  gives  <5C(TS)  =  1 1  GM. 
ea0=  10  a.u.,  £  =  24.55,  and  n2=1.847  (EtOH). 

Experimental  TPA  cross  section  measured  in  EtOH  (relative  error  of  about  1/2  of  magnitude)  (Ref.  16). 
ga0=ll  a.u.,  £=24.55,  and  n2=  1.847  (EtOH). 

hA/z,ge=2.9  D,  using  A/rge=7.3  7.8  I)  (Ref.  62)  (electro-optic)  gives  i?g(TS)=50-57  GM;  and  using  Ayu,ge 
=  2.72  D  (Ref.  15)  (solvatochromism)  gives  <5C(TS)=7  GM. 

‘a0=ll  a.u.,  £=32.63,  and  n~- 1.758  (MeOH). 

Experimental  TPA  maximum  cross  section  measured  in  MeOH  (Ref.  19). 
kAfj.gs=3.0  D,  using  A/zgc=2. 1  I  D  (Ref.  67)  (solvatochromism)  gives  <5C(TS)=4  GM. 

lA/x:,,=4.4  D,  using  A/i,ge=1.7  D  (Ref.  67)  (solvatochromism)  gives  <5C(TS)  =  2  GM;  using  A/xgc=2.8 1  D  (Ref. 
64)  (solvatochromism)  gives  <5C(TS)  =  6  GM;  and  using  A^tge=4.35— 7.01  D  (Ref.  5)  (electro-optic)  gives 
<5C(TS)  =  14-35  GM. 

“Experimental  TPA  maximum  cross  section  measured  in  MeOH;  r/  is  the  quantum  yield  (Ref.  17).  ry=0.6 1 
(Ref.  58),  0.56  (Ref.  61)  in  EtOH.  t?=0.90  in  50%  EtOH  (Ref.  9). 


strengths  are  slightly  smaller  than  the  computed  values  of 
about  0.4-0. 5,  with  those  predicted  by  CAMB3LYP  and 
mCAMB3LYP  functionals  being  the  larger  ones.  The  transi¬ 
tion  dipoles  and  oscillator  strengths  obtained  with  the  PCM 
are  also  slightly  larger  than  the  corresponding  SCRF-S  val¬ 
ues.  The  overall  errors  for  PCM,  however,  remain  low  for 


B3LYP  (0.41  D)  and  PBEO  (0.55  D).  The  MAEs  using  the 
SCRF-S  model  are  0.28,  0.48,  and  0.71  D  for  the  B3LYP, 
»iCAMB3LYP,  and  CAMBLYP,  respectively.  It  is  encourag¬ 
ing  that  the  predicted  mCAMB3LYP  transition  dipoles  are  in 
accordance  with  those  of  B3LYP  and  PBEO.  The  predicted 
changes  in  /xgc  from  cyclohexane  to  MeOH/EtOH  are  small. 
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This  is  consistent  with  the  lack  of  dependence  on  solvents 
for  absorption  transition  moments  of  C102  and  Cl 53,  as  ob¬ 
served  by  Lewis  and  Maroncelli.11 

Upon  alkylation,  the  computed  /j.gc  and  /  values  of  the 
ICT  band  show  a  slight  increase,  as  predicted  by  the  PCM 
using  the  PBE0  (0.07  and  0.61  D)  and  B3LYP  (0.07  and 
0.64  D)  functionals  for  Cl  [see  an  expanded  version  of  Table 
I  in  the  EPAPS  (Ref.  68)  Supporting  Information  Section]. 
This  is  in  good  agreement  with  the  experimental  shifts  of 
0.08  and  0.74  D  for  the  oscillator  strength  and  transition  mo¬ 
ment,  respectively.  Note  that  these  observed  shifts  are  also 
well  reproduced  using  the  SCRF-S  model  with  the  B3LYP 
(0.08  and  0.68  D)  and  the  mCAMB3LYP  (0.10  and  0.74  D) 
functionals  but  not  as  accurate  with  the  CAMB3LYP  (0.12 
and  0.79  D)  functional.  The  alkylation  of  Cl 20  to  create 
C480  produces  a  similar  increase  (0.09  and  0.95  D)  in  the 
experimental  /  and  fige  values  from  the  absorption  data  of 
Fisher  et  al . 16,19  and  of  Reynolds  and  Drexhage.61  However, 
the  corresponding  computed  values  for  the  syn  and  anti  con¬ 
formations  of  C480  show  smaller  increases  from  that  of 
Cl 20,  which  is  overestimated  by  the  continuum  solvation 
models.  Furthermore,  the  two  isoenergetic  syn  conformations 
are  0.5  and  0.2  kcal/mol  lower  in  energy  than  their  anti 
counterparts  in  the  gas  phase  and  in  MeOH,  respectively. 
Thus,  rapid  syn-anti  conversions  may  occur  at  room  tem¬ 
perature,  leading  to  a  different  mixture  of  conformers  and 
corresponding  intensities  for  the  ICT  band.  Similarly,  the  syn 
conformers  of  the  trifluorine  substituted  Cl  53  are 
0.3  kcal/mol  lower  in  energy  than  their  anti  counterparts  in 
the  gas  phase.  This  is  identical  to  the  syn-anti  relative  energy 
reported  by  Cave  and  Castner.  However,  recent  calcula- 
tions  by  Sulpizi  et  al.  ’  for  C153  using  the  nonhybrid 
BLYP  functional  with  a  plane-wave  basis  set  and  pseudopo¬ 
tentials  revealed  that  the  anti  conformer  has  the  lowest  en¬ 
ergy.  The  other  anti  conformer  was  located  at  0.8  kcal/mol 
above  the  lowest  anti  congener.  The  two  syn  conformers 
were  found  to  be  0.8  and  1.1  kcal/mol  above  the  lowest  anti 
structure.  We  found  that  the  syn  and  anti  conformers  of  C153 
are  isoenergetic  in  EtOH. 

The  shifts  in  transition  moment  for  the  trifluorine  deriva¬ 
tives  follow  the  alkylation  trends  discussed  above,  increasing 
from  C157(4.97  D)<C307(5.52  D),  C152(5.31  D) 

<C152A(5.88  D),  and  053(5.71  D).  The  predicted  transi¬ 
tion  moments  (5.23-5.69  D)  for  051  are  significantly 
larger  than  experimental  value  (4.97  D)  in  EtOH  due  to  the 
importance  of  H  bonding.  For  053,  the  predicted  transition 
values  vary  from  5.84  D  (B3LYP)  to  6.40  D  (CAMB3LYP). 
The  corresponding  experimental  value  of  5.71  D  was  ob¬ 
tained  by  using  the  peak  extinction  coefficient  reported  by 
Reynolds  and  Drexhage61  and  the  FWHM  reported  by 
Gustavsson  et  al.  This  value  is  in  good  agreement  with  the 
transition  dipole  moments  of  5.03-5.71  D  reported  by  Lewis 
and  Maroncelli11  for  053  in  methanol,  as  obtained  by  nu¬ 
merically  integrating  the  spectrum  and  scaling  them  with  dif¬ 
ferent  refractive  index  correction  factors.  Their  nonscaled 
value  of  5.80  D  for  the  emissive  transition  dipole  moment, 
which  was  reported  to  be  the  same  as  the  absorptive  moment, 
is  slightly  (0.1  D)  larger  than  the  value  integrated  with 
Gaussian  line  shape  above.  The  predicted  errors  of 


0.3 -0.4  D  for  B3LYP  and  mCAMB3LYP  are  close  to  the 
estimated  experimental  uncertainty  of  0.3  D  (5%)  for 
C153.11 

B.  Two-photon  absorption 

The  TPA  absorption  cross  sections  are  listed  in  Tables 
II-IV  for  various  basis  sets  and  functionals  along  with  avail¬ 
able  cross  sections.  In  addition,  cross  sections  obtained  with 
experimental  /j.ge  (in  Table  I)  and  the  difference  between  the 
ground  and  first  excited  state  dipole  moments  (A/rge)  using 
the  two-state  approximation  are  also  listed  for  four  cou- 
marins.  These  experimental  A/zgc  have  been  estimated  by 
several  experimental  methods,  including  time-resolved  mi¬ 
crowave  dielectric  absorption  measurements,62  which  are 
claimed  to  be  superior.  Note  that  these  experimental  studies 
were  not  carried  in  the  same  solvents  as  in  the  TPA  studies. 
Thus,  the  cross  sections  obtained  from  these  quantities  offer 
qualitative  rather  than  quantitative  comparisons  with  the  TPA 
measurements.  The  consensus  in  OPA  intensity  line  shape  at 
room  (or  higher)  temperature  appears  to  be  Gaussian.  How¬ 
ever,  both  Gaussian  and  Lorentzian  functions  are  often  used 
to  compute  TPA  cross  sections.  The  difference  between  the 
maxima  of  the  two  line  shape  functions,  as  shown  for  Cl 20 
in  Table  II,  is  a  factor  of  1.48,  with  the  Gaussian  cross  (<5C) 
section  being  larger  than  the  corresponding  Lorentzian  (<5L) 
value.  TPA  cross  sections  computed  with  the  single  residue 
of  the  quadratic  response  function  are  denoted  by  SRQR. 
DRQR-SOS  refers  to  values  computed  with  the  SOS  method 
using  the  DFT  ground  state  dipoles,  transition  dipoles  and 
excitation  energies  from  linear  response  TDDFT,  and  the  ex¬ 
cited  dipoles  and  excited-to-excited  transition  dipoles  from 
the  double  residue  quadratic  response  function.  The  latter 
approach  is  computationally  more  demanding,  but  offers  in¬ 
sights  into  the  mechanism  of  TPA.  For  Cl 20  with  the  lowest 
five  excited  states  in  the  SOS  expression,  the  two  methods 
yield  cross  sections  at  the  TPA  resonance  within  a  few  GM 
of  each  other  (see  Table  II).  Our  SOS  calculations  reveal  a 
type-II  TPA  for  the  first  state  of  Cl 20,  in  which  the  TPA  is 
driven  by  A/rge.  Thus,  the  difference  in  the  cross  sections 
obtained  from  the  two-state  and  the  five-state  approximation 
is  (1-3  GM)  small  (cf.  Tables  II  and  III).  We  found  that  the 
cross  sections  are  not  significantly  affected  by  the  basis  set. 
For  C120,  changes  of  Sc  upon  going  from  the  6-31G(c/)  to 
the  6-311  +G(d,p)  basis  set  using  the  SRQR  method  are  less 
than  1  GM  for  all  functionals.  For  the  largest  basis  set,  varia¬ 
tions  in  the  peak  cross  sections  are  within  2  GM,  as  obtained 
with  CAMB3LYP  (5.5  GM),  PBEO  (6.9  GM),  and  B3LYP 
(7.1  GM)  functionals.  Small  variations  among  the  function¬ 
als  are  also  observed  in  solvents  for  Cl 20  and  other  systems. 

In  solvents,  the  predicted  cross  sections  of  13-14  GM 
for  C120  (Table  IV)  are  significantly  larger  than  those  in  the 
gas  phase,  due  to  the  increase  in  A/zgc  (Table  III).  The  cross 
sections  predicted  by  the  PCM  are  significantly  larger  than 
those  obtained  with  the  SCRF-S  using  the  two-state  approxi¬ 
mation.  This  can  be  attributed  to  larger  excited  state  dipole 
moments  predicted  by  the  PCM.  It  is  not  straightforward  to 
compare  these  computed  values  with  experiment,  since  the 
TPA  cross  sections  of  19  and  14  GM  measured  at  3.29  and 
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4.13  eV  in  ethanol16  and  water,18  respectively,  are  not  the 
local  maxima  of  the  TPA  spectra.  The  maxima  are  likely 
located  at  higher  (lower)  energy  with  even  larger  cross  sec¬ 
tions  than  the  predicted  values  for  ethanol  (water).  Interest¬ 
ingly,  these  cross  sections  fall  in  the  range  (2-26  GM)  of 
values  obtained  for  the  two-state  approximation  using  ex¬ 
perimental  yu,ge  and  A/U,ge. 

The  predicted  TPA  cross  section  (energy)  increases  (de¬ 
creases)  with  an  increasing  degree  of  alkylation.  This  is  con¬ 
sistent  with  increase  in  A/rgc  and  /rgt  upon  alkylation  that  has 
been  predicted-9  and  experimentally  observed.62  However, 
the  increase  in  the  cross  section  upon  an  increasing  degree  of 
alkylation  (C120<C1  <  C480)  is  not  quantitatively  repro¬ 
duced.  The  experimental  cross  sections  of  104  and  160  GM 
for  Cl  and  C480  (Table  IV)  are  significantly  larger  than  the 
corresponding  predicted  values  of  28  (Cl)  and  22 
(C480)  GM.  These  discrepancies  can  be  attributed  to  the  un¬ 
certainty  both  in  the  predicted  and  observed  cross  sections. 
In  addition  to  the  errors  from  the  linewidth  functions  as  dis¬ 
cussed  above,  the  computed  errors  might  be  attributed  to  the 
specific  solvent  interactions  that  are  not  accounted  for  in  the 
SCRF  model.  Experimentally,  the  errors  for  TPA  cross  sec¬ 
tions  for  Cl  and  Cl 20  have  been  reported  of  about  half  order 
of  magnitude.16  For  C480  (C102),  the  TPA  cross  section  was 
obtained  by  relative  fluorescence  measurement19  assuming 
unity  fluorescence  quantum  yield  for  both  the  p-bis- 
o-methyl-styrylbenzene  (BMSB)  reference  and  C480,  which 
was  reported  to  be  0.58  (Ref.  61)-0.95  (Ref.  9)  in  EtOH. 
Very  high  (0.995-0.870)  quantum  yields  were  reported  for 
BMSB  in  a  few  nonpolar  solvents,  but  was  not  given  for 
EtOH.6'  Thus,  an  adjustment  the  TPA  cross  section  of  C480 
based  on  measured  quantum  yields  would  be  useful. 

For  the  trifluorinated  series,  the  predicted  TPA  cross  sec¬ 
tions  follow  the  alkylation  trends  previously  observed.  The 
TPA  maximum  ( SG )  cross  section  (energy)  is  predicted 
(SCRF-S-BSLYP/bd  1G '"")  to  increase  (decrease)  from 
14.0  GM  (3.45  eV)  to  25.1  GM  (3.17  eV)  for  C151  and 
C307,  respectively.  The  TPA  spectrum  of  C151  in  EtOH 
taken  from  2.95  to  3.28  eV  shows  a  cross  section  of  47  GM 
at  3.28  eV,16  which  is  comparable  to  the  OPA  maximum  of 
3.24  eV  in  the  same  solvent.  The  experimental  TPA  cross 
section  of  47  GM  for  C151  is  more  than  twice  maximum 
cross  section  of  19?7±5.5  GM  obtained  in  MeOH  at  3.20  eV 
for  C307,17  even  with  the  fluorescence  quantum  yield  77=1. 
Thus,  the  observed  redshift  in  the  TPA  maximum  of  C151 
upon  alkylation  is  consistent  with  the  theoretical  prediction. 
However,  the  decrease  in  the  corresponding  cross  section  is 
in  discordance  with  theory.  It  should  be  noted  that  the  ex¬ 
perimental  cross  section  of  47  GM  reported  for  Cl 51  has  a 
much  larger  uncertainty  (half  order  of  magnitude)16  than  that 
of  C307.  Thus,  new  experimental  studies  on  the  TPA  cross 
sections  of  C151  and  other  coumarins  would  be  useful  in 
verifying  the  predicted  results. 

IV.  SUMMARY  AND  CONCLUSIONS 

We  found  that  gas-phase  excitation  energies  for  Cl 52, 
C152A,  C522,  and  C153  are  well  reproduced  with  the 
6-3 1 1G(<V,/?)  basis  set,  with  the  mean  absolute  errors  of  0.03, 


0.08,  0.10,  and  0.32  eV,  for  PBE0,  B3LYP,  mCAMB3LYP 
(a=0.19  and  /3=  0.19),  and  CAMB3LYP  (a=0.19  and  /3 
=  0.46),  respectively.  The  corresponding  sets  of  0.16,  0.05, 
0.06,  and  0.29  eV  are  found  for  the  6-31 1  +  +G(d,p),  which 
has  negligible  differences  in  S0-Sl  excitation  energies  and 
oscillator  strengths  from  the  larger  aug-cc-pVTZ  basis  set. 
All  the  four  functionals  predict  J>0-5j  energy  to  shift  to  the 
red  with  an  increasing  degree  of  alkylation.  Trifluorination  at 
the  4-methyl  group  gives  pronounced  redshifts  that  increase 
in  magnitude  with  a  higher  degree  of  alkylation.  These  pre¬ 
dictions  are  found  to  be  in  agreement  with  the  experiment. 

The  6-31 1G (d,p)  excitation  energies  predicted  by  the 
SCRF-S  model  overestimate  experimental  S0-S\  energies  by 
0.05,  0.22,  and  0.42  eV  for  the  B3LYP,  /«CAMB3LYP,  and 
CAMBLYP,  respectively,  in  cyclochexane.  In  MeOH/EtOH, 
the  MAEs  for  excitation  energy  (transition  moment)  are  0.13 
(0.28),  0.27  (0.48),  and  0.48  (0.71)  eV  (D)  (D)  for  the 
B3LYP,  mCAMB3LYP,  and  CAMBLYP,  respectively.  The 
corresponding  MAEs  for  the  PCM  are  0.11  (0.41)  and 
0.21(0.55)  eV  (D)  for  B3LYP  and  PBE0,  respectively.  The 
larger  errors  of  excitation  energy  predicted  in  alcoholic  sol¬ 
vents  are  attributed  to  the  H  bonding,  which  is  not  accounted 
for  in  the  continuum  solvation  models. 

The  first  excited  state  of  coumarins  has  large  TPA  cross 
sections,  which  are  driven  by  the  change  in  dipole  moments 
between  the  ground  and  the  more  polar  first  excited  state.  We 
observed  a  large  enhancement  in  the  TPA  cross  section  in 
polar  solvents  due  to  the  increase  in  excited  state  dipole  mo¬ 
ments  compared  to  the  corresponding  ground  state  values. 
We  found  that  the  variations  in  the  TPA  peak  cross  sections 
are  small  (1-3  GM)  among  the  various  basis  sets  and  func¬ 
tionals,  including  the  ;«CAMB3LYP  functional  which  sig¬ 
nificantly  reduces  the  errors  in  the  TPA  and  OPA  resonances. 
The  predicted  cross  sections  are  largely  consistent  with  ex¬ 
periment  to  within  the  experimental  uncertainties.  Some 
large  deviations  may  be  attributed  to  the  large  uncertainties 
in  the  experimental  cross  sections.  Thus,  new  experimental 
TPA  studies  for  these  coumarins  would  be  useful  in  verifying 
the  predicted  results. 
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